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ABSTRACT 
In  order  to  determine  the  shape  of  the  ionization  curve  in 
nuclear  emulsions  at  values  of  /f>100s  we  have  measured  the  blob 
density  of  relativistic  electron  tracks.   The  pv  of  the  pair  pro- 
duced electrons  used  was  determined  by  multiple  scattering.   A 
total  of  84.45  centimeters  of  electron  track  was  bl#b  taunted   and 
scattered.   From  T  =  100  to  T  -   5400  the  data  were  combined  into 
14  points  each  with  a  statistical  uncertainty  on  blob  density  of 
less  than  1%.   These  points  indicate  a  level  "plateau"  and  show 
no   deviation  from  this  plateau  within  our  statistical  accuracy. 
Pions  were  used  to  estimate  the  minimum  of  the  ionization  curve 5 
and  ratio  of  blob  density  plateau  to  blob  density  minimum  is 
estimated  to  be  1.140  +  .020. 
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1.   Introduction 

The  classical  theory  of  the  energy  loss  of  a  charged  particle 
passing  through  matter  was  developed  by  Bohr  /I/  and  later  ex- 
tended quantum  mechanically  by  Be the  /2/  and  Bloch  /3/.   These 
theories  predicted  that  the  energy  loss  by  ionization  per  unit 

path  length  was  a  function  of  velocity  and  charge.   In  particular 

2 

the  energy  loss  was  to  decrease  approximately  as  1/v  to  a  minimum 

at  v/c  —  0.96.   After  this  minimum  the  energy  loss  should  loga- 
rithimically  increase  indefinitely  due  to  relativistic  effects 

on  the  electric  field  of  the  charged  particle.   Swann/4/  predicted, 
however,  that  the  increase  in  energy  loss  would  be  limited,  in 
condensed  substances,  by  the  polarization  of  the  medium.   Quanti- 
tative calculations  by  Fermi  /5/,  treating  the  electrons  in  the 
medium  as  classic  oscillators  of  one  frequency,  showed  that  the 
most  probable  energy-loss  per  unit  path  length  should  increase 
to  and  remain  at  some  maximum  value  which  depends  on  the  density 
of  the  medium.   This  maximum  value  was  therefore  called  the 
"Fermi  Plateau",  and  the  polarization  effect  the.  "density  effect." 
Fermi's  theory  has  been  extended  and  elaborated  upon,  using 
multi -frequency  theory,  by  Wick  /6/,  Halpern  and  Hall  /7/s 
Schbnberg  /8/,  and  Sternheimer  /9/;  but  these  extensions  show  no 
significant  difference  at  v~  c  and  also  predict  a  "plateau." 
Although  emulsion  groups  were  eager  to  measure  the  relativistic 
and  density  effects  in  solid  media,  it  was  not  until  1948  that 


Berriman  introduced  photographic  emulsions  that  were  sensitive 
to  singly  charged  particles  of  minimum  ionization. 

The  first  attempts  to  measure  the  relativistic  rise  met  with 
negative  results.   These  early  experiments  are  discussed  by  Price 
/10/.   However,  in  1950  Pickup  and  Voyvodic  /ll/  reported  the 
first  experimental  indications  of  a  relativistic  rise.   Corson 
and  Keck  /12/,  McDiarmid  /13/,  and  Morrish  /14/  reported  the 
existence  of  an  ionization  plateau;  but  none  of  these  experiments 
carried  their  measurements  down  to  the  ionization  minimum.   Two 
more  experiments  published  in  1952  showed  a  relativistic  increase 
in  grain  density  from  minimum  to  plateau,  but  differed  in  the 
magnitude.   Daniels  et  al   /15/  reported  a  rise  of  about  8%, 
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while  Shapiro  and  Stiller  /16/  obtained  an  experimental  ratio  of 

+0.04 
1.12-0.03.   Stiller  and  Shapiro  /17/  duplicated  their  experiment 

the  following  year  using  blob  count  instead  of  grain  count  as  a 
measure  of  ionization  and  found  virtually  the  same  results-- 
a  saturation  of  the  curve  at  If  > 100,  and  a  rise  above  minimum 
of  about  1^%.   Michaelis  and  Violet  /18/  and  Fleming  and  Lord 
/19/  also  attempted  to  measure  the  ratio  of  ionization  at  plateau 
to  ionization  at  minimum.   In  both  experiments  twi<  points  were 
plotted:   one  near  minimum  and  one  near  plateau.   These  experi- 
ments have,  underlying  their  design,  the  implicit  assumption  of 
a  constant  Fermi  plateau.   Experimenters  were  now  trying  to 
measure  two  parameters:   at  what  value  of  'Y*  the  "plateau" 


begins,  and  how  far  above  the  ionization  minimum  the  "plateau" 

lies.   The  results  of  Alexander  and  Johnston  /20/  are  quoted  as 

a  value  of  grain  density  plateau/grain  density  minimum  of 

1.133  +  .008,  and  Jongejans  /21/  has  stated  in  1960: 

It  seems  that  a  remarkable  agreement  exists  about 
the  ratio  plateau  to  minimum;  we  find 

g    =  1.129  +   .010 
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This  value  compared  with  1.133  +   0.008  of  Alexander 
and  1.14  +   0.03  of  Stiller. 

Recently  an  extensive  study  of  the  relativistic  rise  of  grain 

density  was  carried  out  by  Patrick  and  Barkas  /22/s  but  their  main 

attention  was  focused  en  the  rise  from  minimum  and  transition  to 

plateau. 

Theoretical  studies  /23-26/  were  continuing,  but  the  major 

effort  was  devoted  to  finding  the  exact  shape  of  the  curve  in  the 

interval  during  the  relativistic  rise  and  transition  to  plateau. 

If  the  ionization  curve  could  be  plotted  accurately  in  this  range 

of  T  ,  the  identification  of  particles  with  velocities  in  this 

range  could  be  made.   The  rate  of  energy  loss  per  unit  path  could 

be  measured;  and  together  with  a  measurement  of  velocity,  would 

give  the  particle  identity  even  on  tracks  which  passed  from  the 

emulsions.   It  is  not  surprising  that  the  plateau  region  was  given 

relatively  little  attention.  Price  /10/  perhaps  best  sums  up  the 

prevailing  feeling  when  he  says: 


This  is  a  field  of  work  in  which  the  theory  was 
largely  established  well  before  any  experimental 
verification  was  possible,  and  in  which  specta- 
cular discoveries  were  neither  expected  nor 
obtained. 

In  1.962,  however,  Alekseyeva  et  al  111  I   reported  such  a 
discovery.   They  announced  that  their  theory  predicts  instead  of 
a  'flat  plateau,  a  decrease  of  the  ionization  loss  at  very  high 
velocities.   Furthermore,  they  stated  that  their  experimental  work 
was  in  agreement  with  their  theory.   Their  data  show  a  drop  of 
several  per  cent  in  the  blob  density  of  electrons  in  the  region 
between  T  -  200  and  T  ^  600  (T—  l/)fi~  V*/t  1      In-   work  done 
at  the  same  time,  the  data  of  Stiller  /28/  show  a  slight  tendency 
to  "peak"  at  T-  750.   Since  previous  data  above   "Y*  SS  200 
had  a  statistical  accuracy  no  better  than  2%,  earlier  experiments 
were  able  neither  to  confirm  nor  to  refute  the  Russian  theory. 
With  these  facts  in  mind,  we  decided  to  attempt  the  task  of  gather- 
ing sufficient  data  to  provide  a  statistical  accuracy  of  1%   at 
high  values  of   *|   .   Since  the  area  of  high   I   had  pre- 
viously been  neglected,  we  decided  to  concentrate  our  efforts  in 
this  area  and  to  take  only  enough  data  below   T   =  100  to 
establish  the  fact  that  our  measurements  were  sensitive  enough  to 
detect  any  variation  in  the  curve.   The  precise  details  of  the 
formulation  of  the  experiment  are  given  in  the  next  section. 

We  have  multiple  scattered  a  total  electron  track  length  of 
84.45  centimeters  and  counted  191>-i95  electron   blobs.   We  find 


no  significant  evidence  for  a  departure  from  a  flat  "plateau"  in 
the  region  I   ^   100. 


2.   Formulation  of  the  Experiment 

After  the  purpose  of  this  experiment  had  been  decided,  it 
was  necessary  to  formulate  our  exact  procedures.   Sine-,  we   wanted 
to  measure  particles  with   '   as  high  as  possible,  we  chose  to 
work  with  electrons.   One  source  of  these  high   »   electrons  is 

pair  production,  using  the  gamma  rays  from  the  decay  of  neutral 
pions  for  initiation  of  the  pair  production. 
-17 

The  neutral  pions  are  obtained  from  interactions  of  high 
energy  primary  beam  particles  with  nuclei  in  the  emulsion. 
K.  I.  Alekseyeva  et  al   /27/  used  this  same  procedure,  using  19.6 
Rev.  protons  as  the  primary  beam  particles. 

The  choice  of  primary  beam  particles  was  made  with  two 
criteria  in  mind.   First,  the  beam  energy  should  be  as  high  as 
possible,  since  the  number  of  neutral  pions  formed  in  the  primary 
star  (and  hence  the  number  of  pair  produced  electrons)  is  a  func- 
tion of  the  beam  energy.   Sec  and,  Shapiro  /29/  indicate.-  that  the 
particles  used  for  normalization  should  have  T  >  100.   To  satisfy 
these  two  conditions,  pellicles  were  used  that  had  an  incoming 
Tf~       beam  of  16.2  +  0.6  Bev.   This  corresponds  to   iT^  115 
and  enabled  us  to  use  the  primary  beam  particles  as  our  normalisa- 
tion particles.   Because  of  the  source  of  electrons,  an  area  one 
centimeter  square  located  seven  centimeters  from  the  primary 


beam  entrance  end  was  chosen.   This  area  was  then  scanned  for 
electron  pairs.   There  were  42  electron  pairs  found  in  pellicle 
fl     -132,  67  in  pellicle  $\    -133,  and  72  in  pellicle  H  -134. 

The  pellicles  we  used  were  from  a  600  micron  Xlford  K-5 
emulsion  stack  which  was  exposed  to  a  16.2  +.6 lev  TT   beam 
at  CERN  in  1960.   This  stack  was  flown  to  Berkeley,  California, 
and  developed  there  four  days  after  exposure.   Of  the  180  pelli- 
cles in  the  stack,  numbers  121  through  180  are  at  the  United 
States  Naval  Postgraduate  School  on  loan  from  Br.  Walter  H.  Barkas 
at  Lawrence  Radiation  Laboratory. 

Although  it  would  have  been  desirable  to  measure  the  rate  of 
ionization  loss  over  the  entire  energy  range  with  a  single  species 
of  particle,  electrons  of  energy  less  than  30  Mev  are  difficult 
to  measure.   This  difficulty  arises  largely  because  at  these  energies, 
electrons  suffer  large  multiple  scatterings.   Because  of  the  strong 
scattering:   (1)  it  is  difficult  to  contain  the  electron  in  an 
emulsion  for  any  significant  length;  (2)  the  probability  of  in- 
advertantly shifting  to  another  track  is  high;  (3)  it  is  difficult 
to  separate  the  background  count  from  the  true  count;  and  (4)  the 
measurement  of  the  electron  track  length  becomes  somewhat  sub- 
jective. We  therefore  used  pions  and  protons  from  the  primary 
stars  as  our  source  of  data  for   !  <,  40.  Although  it  was  not 
our  intention  to  measure  accurately  the  details  of  the  rise  to 
plateau,  we  wanted  to   be  sure  that  our  measurements  were  sensitive 


enough  to  detect  a  rise. 

In  order  to  calculate  the  pv  values  of  the  electrons,  protons, 
and  pions,  a  program  was  written  for  the  CDC  1604  computer  at  the 
U8NPGS.   A  discussion  of  this  program,  along  with  the  program  flow 
chart  and  printout  are  included  it!  the  appendices.   Three  physi- 
cists did  the  multiple  scattering  on  each  track,  scattering  each 
electron  track  at  least  twice  tc  prevent  misidentification.   For 
fast  electrons  the  determination  of  pv  is  complicated  by  the  radia- 
tive energy  loss  along  the  electron  path.   In  :rder  to  minimise 
errors  resulting  from  this  effect,  no  single  electron  oraok  was 
•scattered  for  more  than  1.5  centimeters,  and  most  were  scattered 
for  only  0.5  centimeters.   Each  computer  program  output  was 
analyzed  for  bremsstrahlung  effects  (see  appendices),  and  If  a 
bremsstrahiung  was  suspected,  the  track  was  discarded.   Furthermore, 
the  same  segment  of  electron  track  used  in  the  measurement  of  iv 
was  used  in  the  measurement  of  blob  density.  Thus  if  some  un= 
detected  bremsstrahlung  did  occur  along  the  electron  track,  the 
blob  density  and  the  value  of   »   would  be  averaged  over  the  same 
range  of  pv. 

As  a  measure  of  the  rate  of  ionization  loss,  we  have  used 
blob  density  rather  than  grain  density.   The  greater  information 
content  of  ether  track  parameters  /22,  26/  is  offset  for  near- 
minimum  tracks  by  the  ease  and  accuracy  with  which  blob  densi 
data  can  be  gathered.   Accuracy  in  this  sense  means  reproduci- 
bility.  Price  /10/  has  stated  that  it  is  possible  t        "   :ain 
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results  by  blob  counting  that  vary  by  no  more  than  0.5^4  between 
observers. 

Although  the  rise  from  minimum  to  "plateau"  might  have  been 
different  had  we  used  grain  density,  the  shape  of  the  "plateau" 
should  not  change.   That  is,  if  a  drop  from  "plateau"  exists  when 
grain  density  is  used  as  a  measure  of  ionization,  it  will  also 
exist  when  blob  density  is  used.   Therefore,  all  data  in  this  paper 
is  presented  in  terms  of  blob  density. 

Two  scanners  were  given  the  task  of  blob  counting  all  tracks 
used.   To  prevent  the  need  for  calibration  between  the  two  scanners, 
each  had  two  pellicles  of  her  own.   All  blob  counting  in  each  pelli- 
cle was  thus  done  by  one  observer.  The  electron  and  pion  tracks 
were  normalized  to  tracks  of  the  primary  pion  beam  at  the  same 
depth  in  the  emulsion.   Periodic  recounting  of  random  tracks  showed 
that  several  months  were  required  to  reach  a  satisfactory  reproduci- 
bility of  1%.   Because  of  this  effect,  much  of  the  data  gathered 
in  the  early  part  of  this  experiment  were  eventually  discarded. 
The  data  presented  in  this  paper  are  from  three  pellicles.  Addi- 
tional data  from  the  fourth  pellicle  will  be  published  at  a  later 
date. 


3.   Measurement  of  particle  velocities 

Multiple  scattering  measurements  were  made  on  each  track  to 
estimate  the  product  of  momentum  and  velocity,  pv.   The  velocity 
and  the  corresponding  value  for   u   were  then  calculated.   Scatter- 
ing measurements  were  performed  on  proton,  el-.  :ron,  and  pion 
tracks  using  a  Koristka  R-4  microscope.   The  calculations  were 
accomplished  on  a  CDC  1604  computer;  a  program  description,  flew 
chart,  and  listing  for  which  are  included  in  the  appendices. 

The  same  basic  computation  was  used  to  estimate  the  value  of 
pv  for  all  particles.   The  expression  /3Q/| 

Kco  (S)3/2 

pv  =         {in  Me.) 

573  D 

was  used  for  calculations.   The  scattering  constant,  K  ,  is  chat 

appropriate  for  a  cut-off  without  replacement  at  four  times  the 

true  mean  second  differeo.ee,  D.   The  ceil  length,  s,  is  in  microns. 

The  scattering  constant  was  evaluated  from  the  expression; 

k £    =675(0.090  +  02 IB Logl0  (Ssty 


where  the  constants  are  those  given  by  Scott  /3l/  and  adjusted 
by  Barkas  /30/.   The  equivalent  cell  length,  s',  corrects  for 
the  dependence  of  K   on  particle  velocity  and  charge.   The  same 

CO 

references  supply  the  expression  for  the  equivalent  cell  length; 

s  -  (0.23  +  0.77  v2/c2)(s") 
for  singly  charged  particles. 
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Noise  was  removed  from  the  scattering  measurements  for  each 
track  by  two  distinct  procedures.   The  methods  used  were;  (1)  cell 
length  variation  /32,30/  using  cell  lengths  of  one,  two,  and  three 
times  the  prime  cell  length,  and  (2)  subtraction  of  a  constant 
noise  appropriate  to  each  observer. 

Using  cell  length  variation  the  ncis^  can  be  estimated  from: 


Noise  squared  (i,j)  =  3, 


In  this  expression  \Qk       CtHCI     \I\      are  the  means  of  the  absolute 
values  of  the  measured  second  differences  using  cell  lenghts  of 
i  and  j  times  s.   The  indices  i  and  j  vary  from  one  to  three; 
however    x   sgs  J   .   The  factor  of  0.14  in  the  exponent  comes 
from  considering  the  scattering  constant  as  a  function  of  cell 
length.   The  noise  was  calculated  using  cell  overlap  combinations 
of  (1,3),  (1,2),  and  (2,3)  times  the  prime  cell  length.   The 
three  values  were  compared  for  consistency  but  the  value  of  .noise 
using  cell  overlap  of  (1,2)  was  the  only  one  used  in  estimating 
th  e  true  s  i  gna 1 . 

Repeated  scattering  measurements  were  performed  on  the  16 
Eev  primary  pions  by  all  observers.   For  picns  of  this  energy, 
apparent  scattering  is  almost  completely  noise.  Using  this 
noise,  q,  the  true  mean  second  difference  for  an  electron  track 
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was  estimated  from: 

D2  =    J!2   -  q2 

where  &        is  the  measured  total  mean  second  difference.  Prime 
cell  lengths  of  100 s  250 s  and  500  microns  were  used  with  the 
resultant  noise  found  to  be  nearly  independent  of  these  cell 
lengths . 

All  tracks  were  scattered  twice,  usually  by  different  observer* 
Tracks  'which  did  not  yield  results  in  statistical  agreement, 
either  between  observers  or  method  of  noise  removal s  were  dis- 
carded or  rescattered. 

For  fast  electrons,  the  determination  of  pv  is  complicated 
by  the  radiative  energy  loss  along  the  path  of  the  electron.  To 
minimize  this,  no  electron  track  was  scattered  more  than  1,5  cm. 
Further,  any  electron  track  which  showed  evidence  of  a  change  in 
pv  along  the  scattered  length,  either  as  a  noticeable  single 
scatter  or  as  indicated  by  the  computer  program,  was  discarded. 
The  computer  divides  each  track  into  segments  ten  cell  lengths 
long  to  facilitate  the  detection  of  changes  in  pv  possibly  caused 
by  bremsstrahlung.   (See  app.  I) 

Of  a  total  of  362  electron  tracks  (181  pairs)  located  in 
the  scanned  areas  of  the  three  plates,  225  were  discarded  because 
of  detectable  bremsstrahlung  or  because  the  length  of  track  in 
one  pellicle  did  not  allow  a  minimum  of  20  prime  cell  lengths. 

In  addition  to  scattering  each  crack  twice,  requiring  agree - 
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merit  between  methods  in  determination  of  D9  and  requiring  that  no 
detectable  bremsstrahlung  occur  in  the  used  portion  of  the  track; 
all  tracks  used  met  the  conditions  shown  in  Fig.  1. 


Mir?. 

Mir?  No, 
of  S 

S  used 
(m'/cron) 

Optimism 
No,  or  S 

p/  4  zoo 

1%  o 

20 

100 

40 

200  <  py  <  1000 

1.5 

20 

ISO 

40 

pV  >  1000 

AS 

20 

soo 

30 

Fig.  1   Useability  Criteria 

Fig.  3  /33/  was  used  to  estimate  the  minimum  cell  length  to  use  for 
scattering  to  expect  a  signal  to  noise  ratio   ^   1.5  for  a 
specific  value  of  pv. 

The  standard  deviation,  0"  ,  of  pv  determinations  was  esti- 

K 


mated  using   0" s /4SfT  /32/j  where  n  is  the  number  of  prime  cells, 
/\    is  the  degree  of  overlap,  and  f  Is  taken  from  Fig.  2. 


?\ 

J 

2 

J 

f 

J 

J.022 

J. 046 

Fig.  2    rerlap  Form  Factors 

Using  the  criteria  in  this  section,  the  estimated  errors  in  pv 
ranged  from  15  to  30  percent. 
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4.   Measurement  of  ionization 

Blab  density  has  been  used  as  the  measure  of  ionization  in 
this  experiment.   The  greater  information  content  of  other  track 
parameters s  such  as  grain  density  or  mean  gap  length  /34/s  is 
offset  for  near  minimum  tracks  by  the  ease  and  accuracy  with  which 
blob  density  data  can  be  gathered.   Bauseh  and  Lomb  mlscroscopes, 
with  magnification  of  about  16  x   100 ,  equipped  with  coordinate 
stages  by  the  Lawrence  Radiation  Laboratory  at  Berkeley,  California, 
were  used  for  blob  counting. 

The.  energy  of  fast  electrons  decreases  exponentially  with 
increasing  track  length  in  the  emulsion.   As  the  blob  density  is 
also  a  function  of  energy,  care  was  taken  to  assure  the  scattering 
measurements  for  second  differences  and  the  average  blob  densities 
were  obtained  over  the  same  portion  and  length  of  each  electron 
track.   Therefore,  if  some  undetected  radiative  energy  loss  occurred 
along  an  electron  track,  the  ionization  loss  and  the  value  of   I 
are  averaged  over  the  same  range  of  pv„ 

Blob  count  calibration  curves  for  the  three  plates  used  are 
shown  in  Fig.  4.   These  curves  describe  the  variation  of  blob- 
density  with  depth  in  the  emulsion  for  each  of  the  pellicles.   To 
obtain  these  plots,  the  pellicle  thickness  was  divided  into  tenths 
and  then  the  blob  density  of  che  primary  pions  within  each  tenth 
of  pellicle  thickness  was  determined.   Each  point  represents 
about  12  tracks  or  5,000  blobs.   These  tracks  were  all  located 
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in  the  same  region  as  the  electron  pairs.   The  variation  in  these 

curves  for  pellicles  processed  together,  clearly  si     .      the  need. 

for  careful  calibration  of  each  pellicle.   (See  subhead  79   Critique] 

All  blob  density  data  were  collected  by  tw<  scanners.   All 
counting  in  a  given  plate  was  done  by  one.   Both  scanner a  were 
required  to  return  periodically  to  their  "normalization  tracks" 
(primary  pions)  to  check  for  any  change  in  subjective  criteria  used 
in  blob  counting.   Changes  did  occur s  and  much  of  the  data  gathered 
in  the  early  part  of  the  experiment  had  to  be  discarded  because 
of  this.   Several  months  were  required  to  reach  reproducibility  of 
about  one  percent.   It  should  be  noted  the  two  scanners  had  a 
continued  systematic  difference  of  about  3%  in  their  counts  of  the 
same  tracks.  As  a  consistency  check  the  shape  of  one  of  the  curves 
of  Fig.  4,  that  for  pJ-134,  was  reproduced  by  each  scanner. 

The  electron  and  pion  blob  densities  were  normalized  by  form- 
ing the  ratio  of  blob  density  of  a  track  to  that  of  the.  primary 
pions  at  the  same  depth  in  the  emulsion.   To  select  the  primary 
pion  blob  density  to  be  used  for  normalizing  a  given  track,  the 
initial  and  final  depths  in  the  emulsion  for  that  portion  of  the 
track  used  in  data  gathering  were  determined.   The  calibration 
curves  were  considered  as  linear  sections  and  a  weighted  average 
track  depth  was  used  to  establish  the  depth  in  the  emulsion  of  the 
primary  pion  to  be  used  for  normalization. 
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FIG.   4 

1  7 


5.   Results 

The  effect  which  this  experiment  was  designed  to  detect  was 
a  5%  variation  in  ionization  loss  at  values  of  "j       greater  than 
200.   For  IT  greater  than  10,  data  were  processed  to  provide  a 
statistical  uncertainty  on  all  points  plotted  in  Fig.  12,  13,  and 
14  of  2%   or  less.   At  lower  rV    values,  no  statistical  uncertainty 
exceeds  3%. 

To  provide  this  degree  of  certainty,  several  tracks  were 
combined  to  produce  each  data  point,  but  each  point  represents 
only  one  type  of  particle.   Tracks  were  added  together,  beginning 
with  the  lowest  value  of  I       ,  until  the  sum  of  the  blob  count, 
N,  for  each  point  was  such  that  100x/ifSJ/-2N  =  2%  or  3%  as  des- 
cribed above. 

Since  the  16.2  BeV  pions  (    ~  115)  to  which  the  stack  was 
exposed  have  an  ionization  loss  which  differs  by  less  than  1% 
from  the  ionization  loss  of  the  plateau,  the  blob  count  for  these 
tracks  was  selected  as  the  value  to  which  all  others  should  be 
normalized.   The  effect  of  normalization  on  the  data  for  each  plate 
is  illustrated  by  comparing  Figs.  6,  7,  and  8  with  Figs.  9,  10, 
and  11  for  plate  numbers  132,  133,  and  134  respectively. 

The  following  relations  describe  how  the  normalized  mean 
blob  count,  c,  the  ordinate  of  Figs.  12,  13,  and  14,  was 
determined  for  each  point. 
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N  =  Number  of  blobs  in  the  ith  track 

1 

N  =  Number  of  blobs  in  the  same  length  of  16.2  BeV  pion  track 
o 

m.  ■  Number  of  blobs  per  100  microns  of  track 

n  =  Number  of  blobs  per  100  microns  of  16.2  BeV  pion  track 
o 


assuming  N.   has 


n.    =  N./m. 
l          11 

N     +        fl 

c.    =  n./ri     =  N./m.n 

or  c.    = 

i 

i           ^ 

1              1       O             110 

lu  .    n 

1      o 

a  statistical  uncertainty  of    yN.  . 


(C    (c.  )  =  N./m.n 
u   N  1        1   1  o 

The  c  is  the  average  of  the  c.'s  weighted  according  to  the  inverse 

square  of  the  statistical  uncertainty.  The  summation  is  over  the 
number  of  tracks  involved.  .  o\ 

,  Ml.  Y^  im 

-  ~    /-  £i^a  2  2 

or   C  —   — — .   Kj         since  (n  /n.  )  adjusted  only  the  third 
decimal  place  of  the  normalized  blob  count.   Therefore  each  track 
is  weighted  according  to  the  number  of  blobs  counted. 

The  value  of    |    for  a  given  point  represents  the  arith- 
metic average  of  the  values  found  for  each  track  represented  by 
that  point.   On  the  plateau,  the  maximum  spread  of  the  values 
represented  by  any  given  point  is  22%  of  the  value  at  that  point. 
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(Fig.  14,  plate  134,  °C       =  2600.)  A  typical  spread  is  about  10%. 
Thus  the  points  are  separated  by  a  distance  somewhat  greater  than 
the  error  in    HP  . 

These  processes  were  used  for  determining  the    [   and  normal- 
ized blob  count  of  each  point  for  each  of  the  three  plates,  sepa- 
rately, to  produce  the  curves  of  Figs.  12,  13,  and  14.   Combining 
these  plates  by  superimposing  the  plateau  regions  produced  the 
curve  of  Fig.  5.   This  method  of  combination  was  dictated  by  the 
overwhelming  statistical  weight  associated  with  the  information 
in  the  plateau  region  as  compared  to  the  information  in  the  re- 
mainder of  the  curve.   (For  electrons,  a  total  track  length  of  about 
84  cm.  was  examined.) 
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6.   Comparison  with  theory  and  other  experiments 

As  was  described  in  Subhead  1,  the  relativistic  rise  to 
plateau  has  been  investigated  by  many  persons.   In  1950,  Picl 
and  Voyvodic  /ll/  reported  indications  of  a  relativistic  increase 
in  grain  density  of  about  107o  to  a  plateau  starting  at    F  «  20, 
Similar  results  were  obtained  by  Daniel  et  al   /IS/.   In  .9  52 
Stiller  and  Shapiro  /17/  reported  a  14  +  3%  rise  to  plateau.   In 
our  experiment,  as  in  Stiller  and  Shapiro's,  blub  densities  were 
substituted  for  grain  densities  to  simplify  counting  and  facilitate 
reproducibility  among  observers.   The  magnitude  of  the  relativi -  . . 
rise  of  blob  count  determined  was  14  +  2%.   in  our  experiment. 

The  rate  of  rise  to  the  plateau  and  the  constant  vales  on 
plateau  indicated  by  our  results  is  essentially  the  same  as  that 
determined  by  Stiller  and  Shapiro  and  many  others  but  extends  the 
investigated  region  on  the  plateau  to  higher  values  of  If  .   The 
effect  noted  by  Alekseyeva  et  al  111  I   was  not  observed. 

Our  data  are  entirely  consistent  with  the  theory  as  determined 
by  Halpern-Hall  /7/,  and  Sternheimer  /9/s  with  the  plateau  begin- 
ning at  HP  greater  than  100.   This  indicates  a  deviation  from 
the  theories  of  Daniel  et  al.  /15/  and  Morrish  /14/  which  predict 
a  plateau  beginning  about   I  =  20. 


31 


7.    Critique 

Although  the  data  seem  to  be  in  good  agreement  with  currently 
accepted  theory,  there  are  possible  sources  of  error  which  could 
obscure  a  departure  from  a  flat  plateau  of  ionization  loss.  These 
are  the  determination  of   i   from  multiple  scattering  measure- 
ments and  the  variation  of  blob  density  with  depth  in  the  emulsion. 

To  check  for  systematic  errors  in  our  scattering  data9  the 
pellicles  were  taken  to  the  Lawrence  Radiation  Laboratory  at 
Berkeley  and  a  random  sample  of  tracks  was  scattered  on  the  Koristki 
MS-2  belonging  to  the  Barkas  group.   Consistent  agreement  within 
statistical  error  was  found,  and  we  consider  this  the  best  evidence 
we  have  for  thinking  that  our  estimates  of  X  are  not  systemati- 
cally in  error. 

Since  our  data  were  insufficient  to  allow  restricting  all 
tracks  to  a  small  fraction  of  the  emulsion  thickness «,  a  correction 
was  required  to  compensate  for  the  variation  of  blob  density 
with  depth.   (See  subhead  4.)  Although  pains  were  taken  in  the 
normalization,  one  must  recognize  that  the  variation  with  depth 
is  a  severe  handicap.   Some  encouragement  comes,  however,  from 
the  consistent  results  obtained  from  pellicles  with  rather 
different  normalizing  curves. 
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APPENDIX  I 


SCATTERING  PROGRAM 


This  program  was  written  for  the  purpose  of  calculating  the 
pv  of  a  singly  charged  particle  in  a  nuclear  emulsion.   The  method 
of  multiple  scattering  is  described  in  the  section  on  determina- 
tion of  velocity.   Our  program  was  written  in  Fortran  60  for 
specific  use  on  a  Control  Data  Corporation  1604  computer.   A  dis- 
advantage for  use  as  a  general  scattering  program,  is  that  in.  order 
to  compute  an  equivalent  cell  length,  an  estimate  of  v/c  must  be 
made.   The  data  printout  format  was  designed  so  that  the  actual 
computer  paper  could  be  filed  and  kept  as  a  record.   During  the 
course  of  the  experiment,  printouts  were  added  when  deemed  advan- 
tageous; until  in  its  final  form,  the  program  outputs  sufficient 
data  to  analyze  the  run. 

The  following  are  the  program  inputs.   On  the  first  card   E 
each  run: 

Columns  1-3  Pellicle  number 

Columns  5-7  Event  number 

Columns  9-10  Prong  number 
*Column  11    Scanner  identification 

Columns  12-17  Date  of  scattering 

Columns  18-22  Primary  cell  length 

Columns  23-25  Number  of  Y  coordinate  readings 
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Columns  26-30  Estimate  of  v/c 
Columns  31-35  Microscope  calibration 

*  This  number  is  used  for  filing  and  for  choosing  a  personal 
noise  reading  written  into  the  program. 

On  the  second  card  of  each  run  are  the  following: 

Column  1   Number  one  if  final  track  for  this  computer 

run;  otherwise  number  zero. 
Column  2   Number  one  if  computer  is  to  calculate  pv 
with  both  calculated  and  personal  noise 5 
otherwise  number  zero. 
Column  3   Always  zero. 
On  the  following  data  cards  there  are  eight  coordinate  readings  to 
each  card,  one  reading  every  ten  columns. 

After  reading  the  input  cards,  the  computer  calculates  a  noise 
reading;  and  using  this  noise  reading,  a  value  of  pv.  After  a 
first  run  the  program  selects  the  personal  noise  reading  for  the 
individual  doing  the  scattering,  returns  to  the  beginning,  and  re- 
calculates pv  using  this  personal  noise  reading.   The  personal 
scanner  noise  was  determined  by  having  each  person  scatter  the 
16.2  Bev  pions  of  the  primary  beam.   Cell  lengths  of  100,  250,  and 
500  microns  were  used  to  scatter  the  pions;  and  since  at  these 
cell  lengths  the  signal  from  a  16.2  Bev  particle  is  virtually  non- 
existent, any  reading  is  noise.  Enough  tracks  were  scattered  by 
each  person  to  determine  a  mean  value  of  individual  noise.   The 
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noise  determined  by  this  method  was  found  to  be  independent  of  the 
cell  length  used. 

The  program  outputs  were  selected  so  that  crosschecks  could 
be  made.   The  outputs  were: 

1.  The  number  of  Y  coordinate  readings. 

2.  The  Y  coordinate  readings. 

3.  All  calculated  second  and  third  dif ferences--if  a 
second  or  third  difference  was  cast  out  by  the  guillo- 
tine method,  the  corresponding  second  or  third  dif- 
ference will  print  out  as  zero. 

4.  Number  of  readings  cast  out  by  the  guillotine  method-- 
This  output,  together  with  output  3,  can  sometimes 
pinpoint  a  bremsstrahlung  along  an  electron  track. 

5.  Mean  second  and  third  dif ferences--The  calculation  of 
second  and  third  differences  is  made  for  primary  cell 
lengths,  double  overlapping  cell  lengths,  and  triple 
overlapping  cell  lengths.   There  are  six  corresponding 
values  of  mean  differences  and  cast  outs. 

6.  The  three  calculated  values  of  noise --Noise  was  calcu- 
lated using  primary  and  double  cell  lengths,  primary 
and  triple  cell  lengths,  and  double  and  triple  cell 
lengths.   Only  the  first  value  was  used  to  compute  pv9 
but  the  other  two  were  printed  out  as  a  consistency 
check. 
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7.  The  actual  second  difference  signal --This  value  is  the 
square  root  of  the  mean  second  difference  squared  using 
primary  cell  lengths  minus  the  calculated  noise  squared 
using  primary  and  double  cell  lengths. 

8.  The  signal  using  the  mean  third  difference  and  the 
ratio  of  outputs  7  and  8. 

9.  The  three  values  of  pv  (using  the  three  mean  second 
differences),  the  estimated  statistical  error,  and  the 
signal  to  noise  ratio,  all  using  calculated  noise. 

10.  The  three  values  of  pv,  the  estimated  statistical 
error,  and  the  signal  to  noise  ratio,  all  using  personal 
scanner  noise--If  outputs  9  and  10  did  not  agree  within 
the  estimated  statistical  error,  the  track  was  either 
rescattered  or  discarded. 

11.  The  mean  second  differences  for  primary  cell  lengths 
divided  into  segments  of  ten  consecutive  second  dif- 
ferences --This  output  was  found  to  give  the  best 
indication  of  a  bremsstrahlung  along  an  electron 
track.   A  bremsstrahlung  can  be  suspected  when  the 
mean  second  difference  increases  to  a  higher  value 
between  successive  segments,  and  then  remains  at 
this  higher  value  in  succeeding  segments. 

There  is  one  noteworthy  aspect  to  this  program.   In  some 
situations  the  calculated  noise  can  be  a  negative  number.   Since 
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this  situation  has  no  physical  significance,  the  calculated  noise 
is  set  equal  to  zero.  When  the  computer  is  later  asked  for  a 
ratio  of  signal  to  noise,  it  can  be  troublesome.   This  situation 
can  be  programmed  out,  or  it  can  be  handled  as  we  did.   When  the 
calculated  noise  from  primary  and  double  cell  lengths  prints  out 
as  negative,  the  signal  to  calculated  noise  ratio  is  ignored-- 
it  is  always  infinite. 
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APPENDIX   II 


A 


OGRAM   JCAT 


"Read 

TlATE       No.  OF ?OlMT5 

Event    Cell  Length 
Scat       Beta 
"Date     calib 


~f?EAD 

Last    Yfl 

^^         I*  I  M 

Ite^-p     lW 


II -0 
Nl-N 

Delta =0. 


LOW  LHART 

O— 


A 


Nstop=N-Z*M 


L-I+Z#M 


N-NI 


Nl-NH 

ri-in-i 

DELTA=Y(r)+Y(l-|) 


"DTFFZKD- 


"Delta 


v 


A 
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TOTAL-Q 
NIX  =  0 

KNIXM  =  0 


"Dear(n\VTotal//\ 


DiFr^lVO 
NIX-NIXtl 

Knik(mV-Knix(4i 


3? 


V 


A 


YE5 


NST0P=N-3*M 


DiFF3(M,r> 

Yfcti*HK^V(DM) 

-Yffi 


-O- 


A 


-o 
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Total =0. 
N1IX  =  0 

KNIX(M+3)  =  0 
NST0P  =  N-3*M 


B-NSTOP-KNIXM 
DBAR3(M>ToTAiyk 


NIA--0 
TOTAL=0. 


"DIFF3(MA-0 

NIX-NIX+I 

KNIX(M+3)* 
KNIXJM43H1 
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A=^i/logf(5.*s) 


A 


A 


FN-L 
FH=M 


A 


Ekr5q(m,lVDearz(m^x2.  - 
(  (dear2(0**2.-  jy&miHmi. 
\    (FN/FM^(AJ-I. 


Rerr5q(m,lV 
(Errsq(m,i3)* 

(CAUB**Z.) 


£> 


-o 
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KHEK(l)-0 


5TG5Q(l)  = 
"DbARlfflxxZ.- 
Errsq(U) 


E55=5 


Es5=Z.*5 


-o- 


Next  \ 


Tace 


/ 


-o 


-t>- 


ESS  I  =  ESS/ 
(.25+.77*E>eta**z) 


FCO(l)=SQRTF(b75*  (.09  + 
.£7£*Lo&IOF(5jfES5l))) 


Ess-a*s 


RATIO  (l> 

SQftTF(<5I&SQ(f)/ 


5IGMICffi  = 

CALIBx- 

5QRTF(5I&!>Q(I)) 


A 


■o- 


PBCCO=(FCO(I> 

E^**l.5)/(573.* 
SIGMIC(I)) 


KHEK(I)--+I 


-o 


A 
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KHEK3-0 


COMPUTE 

5IG3 

KHEK3=+I 


J 


A 


Identify 
Scanner  And 

Personal  Noise 


Print 
Filing 
BATA 


Print 

DlFFZKl) 

DlFF3(M)T) 

CK5TOUT5 


-4^ 


■£>■ 


Print 
TJBAR(I)  I-lToG 
DBARSQ(l) 
NOISE 
JMOISE  "SQUARE^ 


Print 

PBCClU-lTDb 

SIGNAL  TO 

Notse  Ratio 


CALCULATE 

Dbar  in 
Segments  of 
10  Cell  Lengths 


Print 


Segment  Signals 
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APPENDIX  III 
PROGRAM  LISTING 

PROGRAM  SCAT 

PROGRAM  SCAT         SYMBOL  DEFINITIONS 
IDEV    FIRST  3,  PELL  NBR    2ND  3,  STAR  NBR,  LAST  2,  PR3NG  NRR 
IDSC    SCANNER   NUMBER 
IDATE   MONTH,  DAY,  YEAR   DATA  T^KEN 
S       CELL  LENGTH  IN  MICRCNS 
N       NUMBER  GF  READINGS  IN  EV.NT 
BETA    ESTIMATED  BETA  Al  INPUT 
LAST    NON-ZERO  FOR  LAST  EVENT  IN  -*UN 
Y(I)    ORDINATES  IN  MICROMETER  UNITS 
DIFF2(M,I)  IS   I-TH   2ND  DIFF  F< R  SELL  LENGTH  MS 
OIFF3(M.I)  IS   I-TH   3RD  DIFF  FHP  3ELL  LENGTH  MS 
DBAR2(M)    IS  MEAN  ?NP  DIFF  FOR  CELLS  MS,  WITH  CUT-OTF 
DBAR3(M)    IS  MEAN  3RD  CIFF  FOR  CELLS  MS,  WITHOUT  CUT-OFF 
SIGSQ(N)    IS  SQUARE  OF  NOISE  CORRECTED  SEC  DIFF  FOR  CELL  NS , 

USING  NOISE  FROM  SINGLE  ANO  DOUBLE  CELLS 
ESS1    IS  EFFECTIVE  S  FOR  CALC  OF   KCO 
FCO(N)  IS  KCO   FOR  CELL  NS,  IN  FLT  PT 
KHEK(N)   +1  IF  PPC  CALC  FOR  CELL  NS 

0  IF  NG  BECAUSE  NOISE  GREATER  THAN  SIGNAL 
-1  IF  NG  BECAUSE  SIG  LESS  THAN   2  X  NOISE 
KHEK3   +1  IF  3RD  DIFF  CALC  OK,  0  IF  REJECTED 
1*6  PRINT  U7 
U7  FORMAT  (9H1J. N.DYER  //) 

DIMENSION  Y(200),KNIX(6),DIFF2(3, 193  I ,D IFF  3 (3, 1 98 ) , DBAR2 ( 3 ) , 
1DBAR3(:>)  ,  ERRSO (  3,  3  ) ,  S  I  GSG(  3  ) ,  FCO(  3)  ,  KHEM  3  )  ,RDBAR  {  6  )  ,RE  RRSU  (  3,  3  )  , 
2SIGMIC(3),PBC( 3),RBARS2(3),RATIO( 3) , 
3SEGI20) 

READ   IN   DATA 

1  FORMATUIO,  II,  16,  F5.C,  13,  Fj.3,  F5.3) 

2  FORMAT(8FlC.O) 

3  FORMAT  (11,11,11) 

READ  1,  IDEV,  IDSC,  IDATE,  S,  N,  BETA,  CALIB 
READ  3,  LAST,  INKER.  ITEMP 


READ 


2  ,  (  Y  (  I  ) ,  I  =  1  ,  N ) 


REARRANGE  FOR  POSSIBLE  DUPLICATES 

11=0 

N1=N 

PELTA=0. 

DO  3000   1=1 ,N 

IF  (Yd))  3002,5001,3001 

3001  Y( I )  =  Y(  I  )  +  DELTA 
Y( 1-11 )  =  /( I ) 

GO  TO  iOOO 

3002  Nl=  N1-1 
11=  11+1 

DELTA  =  Y( I )  +  Y( 1-1 ) 
3000  CONTINUE 

N  =  N1 

CALCULATE  SECOND  DIFFERENCES 

00  it   '<*=!,  3 

NSTOP  =  N-2»M 

DO  5   1=1, NSTOP 

L  =  I«-2*M 

K  =  I+M 

DIFF2(M,I)  =  Y(U-2.*Y(K)+Y(  I) 
5  CONTINUE 
h    CONTINUE 


46 


CALCULATE  MEAN  SECOND  DIFFERENCES  WITH 

CUTOFF  AT   U   DBAR 
PO  6   M*l,5 
TOTAL  =  C. 
NIX  =  0 
KNIX(M)  =  0 
NSTOP  =  N-2«M 
102  DO  7  1=  1, NSTOP 

TOTAL  =  TOTAL  +  ADSF ( D I FF2 ( M, I ) ) 

7  CONTINUE 

A  =  NSTOP  -  KNIX(M) 
DBAR2IM)  »  TOTAL/A 
RDBAR{M)rD3AR2(M)*CALIB 
DO  8   1=  1, NSTOP 

IF  (U.»OBAR2(M)-ABSF(DIFF2(M, I) ) )  9,9,8 
9  DIFF2(M,I)  =  0. 
NIX  =  NIX+1 
KNIX(M)  =  KNIX(M)+1 

8  CONTINUE 

100  IF  (NIX)  6,6,101 

101  NIX  =  0 
TOTAL  =  0. 
GO  TO  102 

6  CONTINUE 

CALC   3RD   DIFFS 
DO  11  M=l,3 
NSTOP=  N-3»M 
DO  10  1=  UNSTOP 
DIFF3(M,I)=Y( I+3«M)-3.*Y(  I  +2«M ) *3.» Y ( I +M )- Y ( I ) 

10  CONTINUE 

11  CONTINUE 

CALCULATE  MEAN  THIRD  C  I  FFERENC E S  WITH 
CUTOFF  AT  U   DBAR 

DO  12  M=l,3 

TOTAL=0. 

NIX  =  0 

KNIX  (M*3)=0 

NSTOP=  N-3*M 
17  DO  13  1=1 , NSTOP 

TOTAL=  TOTAL  +  ABS F ( D I FF3 ( M,  I  )  ) 
13  CONT INUL 

B=  NSTOP-  KNIX(M+3) 

DBAR3(M)=  TOTAL/B 

RDPAR(M+3)=DPAR3(M)»CAL IB 

DO  1U  1=1, NSTOP 

IF  {u.»OPAR3(M)-ABSF(DIFF3(M,I) ))  15,15,1U 

15  DIIF?(H,I»=0. 
NI X=NIX+1 
KNiX(M+3)=KNIX(M*3)+l 

1U  CONTINUE 

IF  (NIX)  12,12,16 

16  NIX  =0 
TOTAL  =0. 
GO  TO  17 

12  CONTINUE 

CALCULATE  NOISE  SQUARED 

A  =  3.+l./LOCF(5.*S) 

DO  20   M  =  1,2 

DO  21   L  =  2.3 

IF  (M-L)  2U,21,2U 
2U  FN  =  L 

FM  =  M 

ERRSQ(M,L)  =  DBAR2(M)»*2.-( ( CBAR2 (L  )» »2.-DBAR2( M ) »*2.  ) / ( (  FN/FM ) »• A- 
11.)) 

RERRSQ(MfL)  =  (E.RRSQ<MtL)  )«(CALIB*«2.  ) 
21  CONTINUE 
20  CONTINUE 
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22 
1122 

122 
30 
31 
33 

341 

340 

34 
342 


32 
55 


>000 
200 

201 


202 
60 

61 

63 

64 

66 

67 

69 

70 

72 
730 
1 

74 

750 
1 

76 
77 
78 


DO  5S 
IF(  IT 
ERRSQ 
GO  TO 
ERRSO 
IF  (E 
FRRSG 
IF  (T 
SIGSO 
GO  TO 
ESS  = 
GO  TO 
ESS  = 
GO  TO 
ESS=3 
ESSl 
FCC(  I 
SIGMI 
PBC(  I 
RAT  10 
KHEM 
GO  TO 
KHEM 
CONTI 
IF(  IT 


.41,3 


CA 
US 

1  =  1,5 
EMP)  12 

(1,2)= 

31 
(  M  ,  N  )  I 
RRSQt It 
(1,2)  = 
RAR2(  I  ) 

'IJ     D 

S 
3U2 

2  .*S 
342 

.•S 

=  ESS/I 
)  =SORT 
C(  I  )  = 
)  =  (FC 
(I  )  =  S 
I)  =  +1 

55 
I)  =  0 
NUE 
EMP)  20 


LCULATE  PBC  FRCM 
ING  NOISE  FROM  1 

2,122,1122 
PERSER 


1  AND  2  CELLS 
AND  2  CELLS 


79 


IFIDBA 

KHEK3 

GO  TO 

SIG3  = 

SIG3=S 

KHLK3= 

GO  TO 

ISCAN= 

GO  TO 

ISCAN= 

GO  TO 

ISCAN= 

GO  TO 

ISCAN= 

GO  TO 

ISCAN= 

GO  TO 

ISCAN= 

PERSLR 

GO  TO 

ISCAN= 

PERSER 

GO  TO 

ISCAN= 

PERSER 

PR  INT 

FORMAT 

10X,3H 

PRINT 

FORMAT 

PRINT 

FORMAT 

RH1DIF 

PLACE 

PRINT 

FORMAT 

PRINT 

FORMAT 

PRINT 

FORMAT 

PRINT 

1) 

FORMAT 


R3(2)» 

=  0 
2C2 

SCRTF 
IG3«CA 

+  1 
(60,6 1 
BHJOHTJ 
72 

OHFRL'D 
72 

BHGAIL 
72 

SHSHt  I 
72 

8HHARR 
72 
BHRUSS 

=  l»38. 
72 
8HSWED 

=  654. 
72 
8HDICK 

=  65U 
73 

(  10X, 
N) 
7U.IDE 

(  10X, 
75 

(  15X,5 
3(1), 7 
MATRIX 
76.Y(  1 

(  10X, 
77,Y(2 
( 10X,F 
78,  Y( 
(  10X.F 
79,Y(4 

( 10XF1 


S  NOISE  SQUARED  FROM  CELLS  MS  AND  NS 
2))  30,31,31 

0. 
»»2.  -  ERRSU  1,2)  )  32,32,33 
BAR2( I )»»2.  -  ERRSQ( 1,2 ) 
U0,5U), I 


.23    +    .77*BETA»*2.) 

F(675.*( .09  +  .2  72»LOGlOF ( 5.»ESS1  )  ) ) 

CALie»SQRTF(SIGSQ(I  )  ) 

0( I)*ESS»»1.5)/(573.«SIGMIC(I)  ) 

QRTF( SIGSO ( I)/ERRSO( 1,2 )  ) 


00,2000,2020 

CALCULATE       SIGNAL    F*0M    3RD    DIFFS 
«2-DRAR3( 1)»*2)    200,200,201 


(  (DBAR3(D«»2.-DBAR3(  1)  ••2.)/(1.5»(  (2. »«A  )-!.))) 
LIP 


,63,64,66,67,69,70),     DSC 


LA 
Y 

0745 

E 
U03 

.408 


•  ~r  w  w 

10H    EVENT    NO., 10X.8HSCANNER    ,     10X,6H      DATE, 10X,5H 


V. ISCAN,IDATE,S,N 

I  10,1  OX, A8,  10X,16,1CX,F5.0, 10X, 13    //) 

H    Y( I) ,8X,8H1DIF2( I ) , 7X , 8H2DIF2 ( I ) , 7X, 8H3DI F2 ( I ) , 
X,8H2DIF3( I ),7X,8H3DIF3( I) ) 

IN    HERE 
) 
Fin.n  j 


8X 
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K  =  ,N-9 

PRINT    B5, (Y( I+U) ,UIFF2( 1 , I *3 ) , D I FF2 ( 2 t  I +2 ) i 0IFF2( 3f I *1 ) « 
1 DI FF3( 1 ,  I +2 ) , DI FF  ^( 2,  [  +  1 ) , C IFF3 ( 3tl  ) , I  =  1  ,K  > 
85    FORMAT ( 10XFlO.Of 5XF 10.0 f5XF 10.0 t5XF 10.0,5XF10.0,5XF10.0,5XF10.0) 

K=N-7 
L  =  N-6 
PRINT  62 f <Y( J*3)tOIFF2( lvJ+2),niFF2(2»J+l)tDIFF2(3tJ)t 

U)IFFM1,J+1),DIFF3(2,J),J=K,L) 
62  FORMAT ( 1  OX F 1 0. 0, 5XF 1 0. 0, 5X F  1  0. 0, 5XF  10.0, 5XF1 0.0,5 XF 10.0  ) 
K  =  N-4 

PRINT  65,  Y(K+2),DIFF2(  l.K+1 ) ,0 IFF2 < 2, K ) . 01 FF3 t 1 , K) 
65  FORMAT ( 10XF10.0,5XF10.0,5XF10.0,20XF10.0) 
L=N-3 

PRINT  68,  Y(L+2),DIFF2( 1 ,L+ 1 ) , D  1FF3 (  1,L) 
68  FORMAT! 10XF10.0,5XF10.0,35XF10.0) 

PRINT  71,Y(N) 
71  FORMAT  M0X,F10.0) 
IAWAY=7HCASTCUT 

PRINT  80, I AWAY.KNIX! 1  ) , KNI X ( 2 ) , KN  IX ! 3),KNIX(4) , KN I  X ( 5 ) ,  KN I  X ( 6 ) 
80  FORMAT (  1  OX , A7, 1 4X, 13, 12X , I  3 , 1 2X , I  3,  1 2X, I  3, 1 2X, I  3, 12X, 13  ///) 
ISIG-6HSIGNAL 
IPBC=3HPBC 
IDBAR=UHDBAR 
OPRINT  86,IDBAR,0RAR2( 1) ,DBAR2 I  2 ) , DB AR2 1 3) ,DBAR3 1 1 ),DBAR3!2) , 
10BAR3(^) 
860F0RMAT! 10X,A4,11X,F10.2,5X,F10.2,5X,F10.2,5X,F10.2,5X,F10.2,5X 
1F10.2  /) 

PRINT  35,  RDBAR( 1 ),RDPAR{ 2 ) , RDBAR I  3 ) ,RDBAR ( 4) ,RDBAR(5), 

1RDBAR(6) 
35  FORMAT  (  1 0X , 7HDBARM IC , RX, F  10 . 3,  5X  ,F  1 0. 3 ,  5X ,  F 10. 3»  5X ,F  10  .  3 ,  5X ," 
1F10.3,5X,F10.3  /) 
IN0ISE=8HN0ISE  SO 
IFIRST=5HM1  N2 
IStC  =  5HM1  N5 
ITHIRD=5HMZ  N3 
DO  43  M«lt J 

RBARS2(M)  =  ( (ORAR2(M)#CALIB)«»2.  ) 
U3  CONTINUE 

J0Y=8HDBARSQMI 

PRINT  Ui4,JOY,RBARS2(  1  )  ,  RBARS2  (  2  ) ,  RB  ARS2  (  3 ) 
44  F0RMAT( 10X, A8 , 5X , F 10 . 5, 5X, F 1 0 . 5 , 5X, F 1 0. 5  ///) 
890PRINT  87,INUISE,IFIRST, ERRSQ! 1 , 1 ) ,1  S EC, ERRSQ I  1 , 3 ) , I THIR 0, 

1ERRSQ(2.3) 
87  FORMAT ( lOX,  A8,  5X, A5 , 2X , F10.3, 5X, A5,  2X ,F 10. 3, 5X, A5 ,2 X,F1 0.3  /) 
ICHECK=6HN0ISE 
J0Y2=8HLRRSQMIC 
v   OPRINT  45,J0Y2,  IFIRST.RERRSG!  1,2),  ISECRERRSQ!  1,3),ITHIRD, 
1RERRS0(2,3I 
.  45  FORMAT! 10X,A8,5X,A5, 2X,F10.M,5X, A5,2X,F10. 4, 5X,A5,2X,F1 0.4  ///) 
00  41  M=l,3 
IF  (KHLK(M))  97,98,97 

97  PRINT  5003,PBC(M),RATIO(M) 

3003  FOPMAT!  10X, 5HPBC , 1  OX, F 1 0. 1 , 24H    SIGNAL  TO  NOISE  RATIO,F10.2  // ) 
00  TO  ul 

98  PRINT  42 

42  FORMAT  ( 10X,3HPRC,  10X,2PHN0ISE  IS  GREATER  THAN  SIGNAL  //) 
41  CONTINUE 

LDIF=OHBY  3DIFF 

IF(KHCK3-0)92,92,90 

90  PRINT  91,ISIG   ,LDIF,SIG3 

91  FORMAT! 10X, A6,A8,6X,  F10. 3  //) 
GO  TO  94 

92  IBA0=4HBARF 

PRINT  93, IS IG   ,LDIF,  IBAO 

93  FORMAT!  1CX, A6rA8,6X,A4  // ) 

94  CONTINUE 

PRINT  120,  SIGMIC!  1  ) 

120  FORMAT  { 10X, 1 5HS IGNAL  BY  20IFF, 5X, F 1 0.3  // ) 
BUDGE  =  SIG3  /  SIGMIC!  1  ) 

PRINT  121,  BUDGE 

121  FORMAT  (10X,27HSIGNAL  RATIO  3DIFF  TO  2DIFF.F10.3  //) 
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CALCULATE  01AR  FROM  SEGMENTS  OF 
10   DIFF2M) 
TOTAL=0. 
KA     =0 
MA  =  N-2 
MI=0 

00  110  1=1, MA 
IF  (OIFF  2  ( li  I  > )   1 11,  112, 1 1 1 

111  TOTAL  =  TPTAL  ♦  ABSF ( D I FF2( 1 , I ) > 
MI  =  MI  +1 

IF  (  MI  -  10  )  113.112,  113 
113  IF  (I  -  MA)  110,  112,1 10 

112  IF  (MI)  11U,  1  10,1 1U 
ITU  KA  x  KA  +  1 

R  =  MI 

SEG  (KA)  =  TOTAL/B 
MI  =  0 
TOTAL  =  0. 
110  CONTINUE 

CALC  SIGNAL  IN  MICRONS  FOR  SEGMENTS 
DO  115  1=1, KA 
IF((SEG(I  )  )*»2.-ERRSQ(  1.2)  )  116.116.117 

117  SEG(  I  )  =  (SORTF(  (SEGd  )  )»«2.-ERRSQ(  1,  2)  )  )»  CALIB 
GO  TO  119 

116  SEG  (I  )  =0. 

119  PRINT  1  IB.  SEG( I  ) 

118  FORMAT  ( 10X,22HSFGMENT  SIGNAL  M ICRDNS, 5X, F7.3  /) 
115  CONTINUE 

ESTIMATING  ERROR  IN  PC 
2020  DIMENSION  PBCE^(3) 
PRIM=N 

00  MOO  1  =  1,3 
TEMP  =  KM  <(  I  ) 
FL  =  I 

FACT1  =  1.4.022MFL-1.) 
FACF2  =  PRIM  -  (2.*FL+TEMP) 
OENOM  =  SQRTF(FACT1»FACT2/FL) 
PBCER(I)  =  PBC( I)»1.37/CEN0M 
400  CONTINUE 

IF(  ITEMP)402, 402, 2030 

402  PRINT  403 

403  FORMATC 10X.45HERROR  IN  MEV  CALCULATED  FROM  PG26U  CERN  NOTES/) 
00  4C6  J  =1,3 

PRINT  404,  PBCER(J) 

404  FORMAT( 10X, 17HPBC  PLUS  OR  M INUS 5X ,F 1 0 . 3/ ) 
406  CONTINUE 

ITEMP=1 

IF  (INKER)  18,10,22 

2030  PRINT  2031 

2031  FORMAT  (10X,24HPPC  USING  PERSONAL  NOISE/) 
on  203:  1=1,3 

PRINT  2032,  PBC(I).  PBCER(I),  RATID(I) 

2032  FORMAT  (10X,  F5.C,10X,  13HPLUS  OR  MINUS  5X,  F5.0, 
110X,  21HSIGNAL  TO  NOISE  RATIO  5X,  F5.2/) 

2033  CONTINUE 

18  IF  (LAST)  48,46,48 
48  CONTINUE 
408  END 
END 


SO 
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